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It has been suggested that the metal-insulator transitions in a number of spinel materials with 
partially-filled t^ g d-orbitals can be explained as orbitally-driven Peierls instabilities. Motivated by 
these suggestions, we examine theoretically the possibility of formation of such orbitally-driven states 
within a simplified theoretical model, a two-dimensional checkerboard lattice with two directional 
metal orbitals per atomic site. We include orbital ordering and inter-atom electron-phonon inter- 
actions self-consistently within a semi-classical approximation, and onsite intra- and inter-orbital 
electron-electron interactions at the Hartree-Fock level. We find a stable, orbitally-induced Peierls 
bond-dimerized state for carrier concentration of one electron per atom. The Peierls bond distortion 
pattern continues to be period 2 bond-dimerization even when the charge density in the orbitals 
forming the one-dimensional band is significantly smaller than 1. In contrast, for carrier density 
of half an electron per atom the Peierls instability is absent within one-electron theory as well as 
mean-field theory of electron-electron interactions, even for nearly complete orbital ordering. We 
discuss the implications of our results in relation to complex charge, bond, and orbital-ordering 
found in spinels. 



PACS numbers: 71.30.+h, 71.28.+d, 71.10.Fd 



I. INTRODUCTION 

Transition-metal spinel compounds AB2X4, where X is 
S or 0, have been for many years the subject of intense 
experimental and theoretical activity. The B sublattice of 
the spinel structure forms corner-sharing tetrahedra, giv- 
ing rise to a geometrically frustrated pyrochlore lattice. 
In general, oxides and chalcogenides of transition met- 
als exhibit complex charge- and spin-ordering that are 
often coupled with orbital-orderingJ^ Coupled orbital- 
charge-spin orderings have been investigated widely for 
compounds of the late transition metal ions with active e g 
orbitals within cubic geometry, such as the cuprates and 
manganites. In contrast, such orderings in spinel com- 
pounds are only beginning to be studied 3 . Spinel systems 
are considerably more complicated than the cuprates or 
manganites due to two distinct reasons. First, the B-ions 
in the spinels, which are in an octahedral environment of 
the X anions, often possess partially-filled f 2g ci-orbitals. 
Examples include V 3+ ions in the vanadates of Zn 4 -^, 
Mn£, and CdZ; V 2 ' 5+ in A1V 2 0^; Ti 3 + in MgTi 2 0^; 
Ir 3.5+ in CuIr 2 S^; and Rh 35 + in LiRh 2 0^ii. The 
threefold degeneracy of the t 2g orbitals, along with the 
weaker tendency to Jahn- Teller ( JT) distortions^ in ti g - 
based systems lead to more complicated physics com- 
pared to e g -based systems. Second, the geometric frus- 
tration in the spinel structure, mentioned above, pre- 
cludes simple orderings 3 . 

Very recently, an exotic orbitally-induced Peierls insta- 
bility has been proposed as the mechanism behind the 
metal-insulator transitions in the spinels like CuIr 2 S4 13 . 
Charge-segregation of Ir into formally Ir 3+ and Ir 4+ , 



accompanied by the formation of an octamer of Ir 4+ 
ions with alternate short and long bonds is observed 
in CuIr 2 S4^. This unusual charge-ordering pattern has 
been explained within the context of an effective one- 
dimensionalization and Peierls instability driven by or- 
bital ordering (00). The spinel structure of CuIr 2 S4 con- 
sists of criss-cross chains of Ir-ions with strongest over- 
laps between orbitals of the same kind (i.e., d xy with d xy , 
d yz with d yz ). Within the proposed mechanism, 00 in 
j r 3.5+ j ong i eac j s t completely filled degenerate d xz and 
dy Z orbitals and effective one-dimensional (ID) quarter- 
filled d xy bands. The latter now undergoes a Peierls in- 
stability that is accompanied by period 4 charge-ordering 
• • • Ir 3+ /Ir 3+ /Ir 4+ /Ir 4+ ■ • • and long-intermediate-short- 
intermediate bonds, as seen experimentally 10 . The 00 
and Peierls instability are thought to be coupled, as op- 
posed to independent. Closely related phenomenologics 
have been proposed to explain the metal-insulators tran- 
sitions in MgTi 2 4 i 3 . and LiRh^O^fii. 

While the proposed scenario^ does provide qualita- 
tive explanations, it raises many interesting questions. 
For example, the bond or charge periodicities that result 
from Peierls distortion in real quasi-lD systems depend 
strongly on the bandfilling. Naively then, one might be 
led to believe that unless the 00 is complete and leads 
to integer occupations of the individual i 2ff orbitals, ex- 
actly commensurate distortion periodicities, as suppos- 
edly observed in CuIr 2 S4^ and MgTi 2 04^, are not ex- 
pected. The exactly period 4 distortion in CuIr 2 S4 and 
MgTi 2 04 is therefore a puzzle. Secondly, the role of 
electron-electron (e-e) interactions, which can have im- 
portant consequences on these transitions, is not clear. 
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For example, it has been suggested that the bond distor- 
tion in MgTi 2 04 causes nearest-neighbor Ti 3+ pairs to 
form a spin-singlet state, giving a drop in the magnetic 
susceptibility at low temperature^. Furthermore, e-e in- 
teractions can strongly affect the distortion periodicities 
in ID chains with fillings different from f electron per 
atomic. Addressing these issues requires explicit calcula- 
tions based on Hamiltonians containing all the necessary 
components, that to the best of our knowledge do not 
exist currently. 

We report here the results of explicit calculations based 
on a simplified model system that displays co-operative 
00 and Peierls instability. Our model system is a 
two-dimensional projection of the pyrochlore lattice, a 
checkerboard lattice with two degenerate directional or- 
bitals per atom. As we discuss in the next Section, the 
model captures the effective one-dimensionalization that 
occurs in the spinel lattice. There exist also subtle differ- 
ences between the our model and the theoretical picture 
that has been employed for the spinels^, that we discuss 
in Section HVl We have considered two different electron 
densities, 1 electron per atom and ^ electron per atom. 
We find that a stable Peierls state with ID order occurs 
for an electron density of 1 electron per atom. Impor- 
tantly, even for significant deviations of orbital occupan- 
cies from integer values, the Peierls distortion is purely 
bond dimerization, as would be true in the precisely i- 
filled ID band. We speculate that this is a consequence of 
interactions between the effective ID chains of our model, 
which are not entirely independent. For the case of i 
electron per atom we did not observe an orbitally-driven 
Peierls state even for integer orbital occupancies in the 
absence of e-e interactions. For the same model in the 
limit of infinite e-e interactions, however, a stable bond 
dimerization occurs. This suggests the second main re- 
sult of our paper, namely that e-e interactions are im- 
portant in stabilizing orbitally-induced Peierls states for 
the case of ^-electron per atom within the orbitally de- 
generate checkerboard lattice. We will discuss possible 
implications of this result for the spinel lattice. 

The paper is divided in following subsections: in Sec- 
tion [IT] we introduce our model Hamiltonian; in Sections 
IIII Al and IIIIBI we present numerical results for 1 and | 
electron per atom respectively; and in Section HVl we con- 
clude and further discuss the relationship of our theory 
to spinels and related materials. 



II. THEORETICAL MODEL 

As discussed above, no quantitative calculations ex- 
ist of coexisting orbital and Peierls order in spinel lat- 
tices. This is due to the huge complexity of these systems, 
which include three dimensionality, triply degenerate t2 g 
metal atom orbitals, possible JT distortions and frustra- 
tion, in addition to e-e interactions. We therefore con- 
struct a simpler model based on a checkerboard lattice 
that is easier to handle, but that is expected to bear simi- 




FIG. 1: (color online) (a) 2D checkerboard lattice with 
orbital-ordering pattern given by the -Hoo term in Eq. [3] 
Filled dots denote metal atom positions, with two orbitals 
per atom, xz and yz. Hopping terms are included along x, y, 
x + y, and y~x directions as defined in Eq.(5j-Eq.[8] Hopping 
along the t x and t y directions, indicated by dashed lines, is 
modulated by the inter-atomic e-p coupling, (b) The (ddn) 
overlap of xz, shown in green (light grey), and yz, shown in 
blue (dark grey), along the x- and y-directions. Note that 
along these directions, non-zero hopping matrix elements ex- 
ist only along x- (xz-xz) or y- (yz-yz) direction, (c) The 
(ddn) overlap of xz and yz orbitals along the y — x direction. 
The overlap along x + y is similar apart from a change of sign 
for the inter-orbital terms. 



larities with the spinel problem. Note that although there 
exists a large body of literature on the consequences of 
frustration within the checkerboard lattice with a single 
atomic orbital per sitei£~— , we are unaware of similar cal- 
culations on the present model which deals with doubly 
degenerate orbitals at each site. An important property 
of orbitally-driven Peierls order in these materials is that 
the 00 leads to Peierls order in several different crystal 
lattice directions^. In order to incorporate orbital de- 
generacy, frustration, and the possibility of Peierls order 
in multiple directions we consider the following Hamil- 
tonian for a checkerboard lattice with doubly degenerate 



metal orbitals at each lattice site. 
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The Hamiltonian in Eq. [T] consists of, (i) Hssh that 
contains the kinetic energy and the inter-ion electron- 
phonon (c-p) coupling (Eq. [2]), (ii) an 00 term Hoo 
(Eq. [3J , and (iii) e-e interaction H ee (Eq. 0} that includes 
short-range e-e interactions within each site. We describe 
each of these terms separately below. 

Hssh includes electron hopping between same as well 
as different orbitals. In Eq. [2 d\^ a creates an electron of 
spin a in the orbital 7 of atom i. 7 and 7' correspond 
to d xz and d yz orbitals, which have lobes that are ori- 
ented perpendicular to each other at each site as shown 
in Fig. [TJ The inter- and intra-orbital hopping matrix el- 
ements tf j are based on Slater-Koster parametrization 
of hopping integrals connecting ti g orbitals 24 . a denotes 
a unit vector along the x, y, x + y, or y — x directions. 
Each bond indicated in Fig. [T] connects two orbitals at 
each metal ion site with two orbitals at another site, and 
is hence written as a 2 x 2 matrix: 
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All hopping integrals are in units of the (ddir) matrix 
element (set to -1) involving the two cZ-orbitals of the 
metal atoms and mediated by the p-orbital of the an- 
ion (not shown explicitly in the figure) in between. In 
the above, the small dd5 hoppings have been neglected. 
Note that electron hoppings along the x (y) direction in- 
volve only the d xz (d yz ) orbitals. Thus the 7r-bonding 
among the orbitals, as opposed to the er-bonding in 
spinels, does not preclude the orbitally-induced effective 
one-dimensionalization 25 . 

The inter-ion c-p coupling in Hssh is written in the 
usual Su-Schrieffer-Heeger (SSH) formal. Here a & is the 
e-p coupling constant corresponding to the bond between 
atoms at i and i + a, A^i +a the deviation of this bond 



FIG. 2: (color online) Pattern of the orbital-ordering and 
bond-order modulation for one electron per atom. Squares 
(circles) represent atoms with predominant xz (yz) orbital 
occupation. Bonds alternate in strength along diagonal di- 
rections, with solid (dashed) lines indicating strong (weak) 
bonds. 



from its equilibrium length, and Kssa the corresponding 
spring constant. We include nonzero e-p couplings only 
along the x and y-directions, in keeping with the Peierls 
distortions in CUIX2S4 and MgTi204 involving only or- 
bitals of the same kind, and assume e-p couplings of equal 
strength (a x = a y = a). We take all spring constants 
Kssh = 1- 

The 00 term Hoo breaks the orbital degeneracy be- 
tween the two orbitals at each individual site. ni 1<T = 
d\ ltJ di 1<y is the number of electrons of spin a on orbital 
7 of site i, and = + n.; 7 ^. The coordinate Qi 
couples to the charge density difference between the two 
orbitals on the site, with a corresponding coupling con- 
stant g and spring constant Koo- We fix Koo to the 
value 1. As written in Eq. [31 the relative phase of the 
Qi at each site is unrestricted. Within the model of ref- 
erence [TH the effective one-dimensionalization is a con- 
sequence of 00 alone and in principle will occur even 
in the absence of Peierls distortion. We have used this 
to determine the preferred Qi mode for both 1 and ^ 
electron per site by calculating the orbital occupations 
in several large lattices with open boundary condition 
(OBC). In a large OBC lattice, the charge densities and 
bond orders far from the lattice edges spontaneously as- 
sume the pattern that would occur in the infinite lattice 
for + coupling limit 27 . An alternate approach to deter- 
mine the dominant 00 mode is to calculate in a periodic 
lattice the energies corresponding to each mode for fixed 
distortion amplitude; the dominant mode is simply the 
one with the lowest total energy. In the present case we 
have performed both sets of calculations for both 1 and \ 
electron per site and have determined that the preferred 
00 mode in both cases is the "checkerboard" pattern of 
Fig.[T](a), which can be parametrized in terms of a single 
amplitude \Q\ with Qi — (— where i x and i y 
are the x and y coordinates of the ith atom. 

The third term in the Hamiltonian, H ee (Eq. [5]), in- 
cludes short-ranged Coulomb repulsions. U (U 1 ) is the 
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on-atom Coulomb repulsion for electrons in same (differ- 
ent) orbitals. While exact diagonalization has been used 
successfully for many ID, quasi-lD, and two dimensional 
(2D) lattices involving both e-e and e-p interactions, in 
the present model with two orbitals per metal site the 
Hilbcrt space is too large to treat any meaningful size 
cluster within exact diagonalization. We will therefore 
consider first the non- interacting (U = U' = 0) sys- 
tem, and then consider the effect of U and U' within the 
unrestricted Hartree-Fock (UHF) approximation, with 
no assumption of the periodicity of the UHF wavefunc- 
tion. Specifically, we replace the interaction terms in H ee 
(Eq.Uby 

U ^2 n iia(n llrJ >) - U y^(^ 7 4.)(^ 7 t) 
U' 

+ — [ («i7<j((^7'<j) + {n iY <j>)) 

- (4 Ya d na + H.c.){4 Ya d na + H.c.) 

- (n i7CT }(ni7' CT ) ~ (nna)(n t ya>) + + H.c.) 2 } 



{rii-fu) and (d\ , a di ia + H.c.) are obtained using a com- 
bination of self-consistency and simulated annealing for 
finding their ground state values. The UHF approxima- 
tion often gives unphysical results for large interaction 
strengths and we will primarily focus on small U and U' . 

We treat the 00 and e-p interactions using a standard 
self-consistent approach derived from the equations 



0(H) 
dQ 



d(H) 
OA 



(9) 



The self-consistency equations derived from Eq. |9] are 
used iteratively given an initial starting distortion. In 
the infinite system the 00 or the bond distortion would 
occur for infinitesimally small coupling constants g and 
a. In finite-size clusters, however, due to the finite-size 
gaps between successive energy levels, nonzero coupling 
constants are required before the symmetry-broken state 
appears. In the following we consider g and a close to the 
minimum values needed for the broken-symmetry state 
to occur. 

We performed calculations for lattices up to 16x16 
(256 atoms with 512 orbitals). The primary difficulty 
in solving Eq. [S] is that because of the large number of 
quantities (Q, A iji+a , and UHF average charge densities) 
to be determined self-consistently, the calculations often 
became trapped in local minima before reaching the true 
ground state. In all cases we have taken care that the 
true ground state was reached. Below we summarize our 
numerical results. These are divided into two subsec- 
tions that discuss average charge densities of 1 electron 
per atom and | electrons per atom, respectively. 
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FIG. 3: (color online) Co-operative orbital-ordering and 
Peierls bond alternation for one electron per atom, with 
U = U' = 0. Calculations are for a 16x16 periodic lattice, (a) 
Majority charge density (see text) in the orbitals forming the 
quasi- ID chains following orbital ordering. Here and in (b) 
circles, diamonds and squares correspond to g = 0.6, 0.8 and 
1.0, respectively, (b) AB along the diagonal chain directions 
(see text and Fig. [5]) as a function of a. (c) Bond alternation 
plotted as a function of charge density. Circles show the effect 
of increasing g with constant a=1.3. Diamonds show effect 
of increasing a with constant gr=0.8. Lines are guides to the 
eye. 



III. RESULTS 

A. 1 electron per atom 

The OO term in our model makes the orbital occu- 
pancy of the xz and yz orbitals unequal at each site. To 
measure the degree of 00 quantitatively, we calculate 
the majority charge density (n + ), defined as the charge 
density in the xz orbitals at +\Q\ sites (see Fig. [Ha)). 
These orbitals form quasi-lD chains in the x direction 
(sites denoted by squares in Fig. [5]) . With one electron 
per site, (n + ) ranges from 0.5 to 1, with (n + ) = 1 indicat- 
ing complete 00, and (n + ) = 0.5 implying the absence 
of 00. Due to the symmetries present in our model, the 
quasi- ID chains along x and y directions are identical-the 
charge density in the yz orbitals at —\Q\ sites (denoted 
by circles in Fig. ^ is also (n + ). 

As an order parameter for the Peierls distortion, we 
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measure the modulation of the bond order, 

£i,i+a, 7 - ^(4+a, CT ,A-,7 + H.C.). (10) 

The bond order we are interested in (Eg . [TO]) is the expec- 
tation value of charge-transfer between orbitals of same 
symmetry belonging to neighboring atoms. The charge- 
transfer is directly coupled to the bond distortion Ai^ +a 
in Eq. [2] and hence for nonzero Aj^+a the charge trans- 
fer across consecutive bonds shows periodic modulation. 
The extent of modulation of B iti+gi ^ is therefore a di- 
rect measure of the SSH distortion strength. As shown 
in Fig. [2] we find bond order modulation along x (y) di- 
rection to involve xz (yz) orbitals only. The modulation 
is purely period-2 (dimerization) with alternating strong 
and weak bonds, and hence we use AB, the difference 
between the calculated strong and weak bond orders in- 
volving orbitals of a particular symmetry, as the order 
parameter for the SSH distortion. Because of symmetry, 
the amplitudes of the bond order modulations involving 
the xz orbitals along the ir-direction, and the yz orbitals 
along the y-directions are identical. 

As discussed above, whether or not a co-operative 
orbitally-induced Peierls instability occurs for (n + ) < 1, 
as well as the periodicity of the resultant bond order 
wave are both important issues. We first consider the 
non- interacting limit (U = U' = 0). The co-operative 
nature of the 00 and bond dimerization is shown in 
Fig. [31 where we show the results of our self-consistent 
calculations for 16x16 lattices. As seen in Figs.^a) and 
(b), the orbitally-induced Peierls state appears only for 
(ri+) > 0.8, with the bond distortion a pure bond dimer- 
ization regardless of (n + ). As expected for a co-operative 
transition, (n + ) increases with g, as seen in Fig.^a). For 
each g there exists an a c beyond which there occur si- 
multaneous jumps in (n + ) and AB (the jump in (n + ) 
becomes progressively smaller as g increases.) The mag- 
nitude of a c decreases with increasing g (see Fig. 02[b)) 
To further show the cooperative effect, in Fig. [2^c) we 
show the effects of (i) increasing a at constant g, and (ii) 
increasing g at constant a. While it is to be expected 
that the orbital order parameter (n + ) increases with g, 
or that the bond alternation parameter increases with a, 
we find that either of the coupling constants enhances 
both and AB. 

Next we consider the correlated case with nonzero U 
and U'. As for U = U' = 0, the orbitally-driven Peierls 
state is again bond-dimerized. In Fig. we plot (n + ) 
normalized by its value in the uncorrelated system as a 
function of U for several values of V. Within UHF, the 
combined effect of U and U' can be to either weaken 
or strengthen the distortion: for fixed V, U tends to 
weaken the OO and the bond distortion, while for fixed 
U, U' strengthens both order parameters. Within the 
UHF approximation for one electron per atom, the effects 
of U and U' cancel exactly when U' — \U . 

From the Hamiltonian, the consequence of U' is to 
minimize the intrasite inter-orbital Coulomb repulsion, 
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FIG. 4: Majority charge density as a function of U and U' 
for one electron per atom, normalized with respect to the 
same quantity for the uncorrelated system. Results shown 
are for for 16x16 periodic lattices with g = 0.8, a = 1.2, and 
Koo = Kssh = 1. Circles are for U' = 0.3(7, diamonds are 
for U' — 0.517, and triangles are for U' = 0.7U. Lines are 
guides to the eye. 

which is achieved by orbital ordering. It is thus not sur- 
prising that U' has the same effect as g in Fig. 2J The 
effect of U, as seen in Fig. [H is however an artifact of 
the UHF approximation. In the case of the strictly ID 
i-filled band chain with one orbital per site, exact diag- 
onalization and quantum Monte Carlo calculations have 
shown that the Peierls bond-alternation is enhanced by 
the Hubbard U— . In contrast, the UHF approximation 
predicts incorrectly that U destroys the bond alternation 
in the above case^. Had we been able to perform ex- 
act diagonalization in the present case, we would have 
found similar enhancement of the bond dimerization by 
U. This would have had a profound effect on our overall 
result, reducing significantly the a c or the threshold 
at which the bond dimerization appears. 

The most important conclusion that follows from the 
above is that an orbitally-induced Peierls instability can 
occur even for incomplete OO ({n + } ~ 0.8), and as 
long as the instability occurs at all, the bond order 
wave is period 2 for (n + ) significantly less than 1. In- 
deed, it is conceivable that the threshold value of (n + ) 
at which the bond dimerization appears can be even 
smaller than 0.8 for nonzero e-e interactions. We have 
found no other periodicity or evidence for soliton for- 
mation in our calculations. An interesting aspect of the 
OO driven bond distortion here is the phase relationship 
between the bond order wave states involving the d xz 
and d yz orbitals in Fig. [2J The short bonds along the 
x and y-directions occupy the same plaquettes, yielding 
a structure that is reminiscent of (but different from) 
the valence bond crystal obtained within the Heisenberg 
spin-Hamiltonian for the checkerboard lattice^. Fur- 
thermore, the bond dimerizations along any one direc- 
tion but on different diagonals of the checkerboard lat- 
tice are strictly "in-phase" . Both of these indicate that 
while the bond dimerizations are consequences of effec- 
tive one-dimensionalization, there exist strong 2D in- 
teractions in between both the criss-cross and parallel 
chains. We ascribe the persistence of the bond dimeriza- 
tion for (n + ) < 1 to the commensurability effect arising 
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FIG. 5: (color online) Co-operative orbital-ordering and 
Peierls bond-alternation for 0.5 spinless fermions per atom. 
Calculations are for a 16x16 periodic lattice, (a) Majority 
charge density in the orbitals forming the ID chains. Here 
and in (b) circles, diamonds and squares correspond to g — 
0.6, 0.8 and 1.0, respectively, (b) AB along the diagonal chain 
directions (see Fig. [2} as a function of a. (c) Bond alterna- 
tion versus majority charge density. Circles show the effect 
of increasing g with constant a=1.3. Diamonds show effect 
of increasing a with constant <?=0.8. Lines are guides to the 
eye. 



from the 2D interactions. 



tained similar results after including U and/or U' within 
the UHF approximation. We conclude that for electron 
density away from 1 electron per atom, orbitally-driven 
Peierls ordering does not occur within our model in the 
non-interacting limit or within the mean-field approach 
to e-e interactions. 

We ascribe the absence of bond-charge distortion here 
to the important role played by the interaction among the 
crisscrossing chains within the checkerboard lattice. We 
have already pointed out in Section MIA that the stabi- 
lization of the perfect period 2 distortion for the case of 1 
electron per atom, even in the absence of complete OO, is 
a signature of such a 2D interaction. Similar 2D interac- 
tions should be relevant also for carrier density 0.5. Since 
within mean field theory the Peierls instability in 2D is 
limited to carrier density of 1, the absence of the bond- 
charge distortion in the present case is to be anticipated. 
On the other hand, we have recently shown in a series 
of papers that specifically for this carrier density and 
one orbital per site, nonzero e-e interactions can strongly 
stabilize bond-charge ordered states in 2D lattice o 27 ' 29 . 
In the case of the checkerboard lattice with a single or- 
bital per site, plaquette spin-singlet formation (though 
without charge-ordering) has similarly been found for the 
same carrier concentration 30 . It is conceivable that sim- 
ilar effects of e-e interactions persist in the present case 
with two degenerate orbitals per site. 

Unfortunately, performing a realistic calculation with 
finite U that goes beyond the UHF approximation in the 
present case is beyond our computational capability. We 
have therefore investigated our model Hamiltonian Eq. Q] 
in the limit of U — > oo, where we assume band orbital 
occupancy corresponding to that for spinless fermions. 
Fig. [5] shows the same order parameters as in Fig. [3] for 
the spinless fermion case. While the transition occurs 
here for a slightly larger value of the coupling constant 
a, it is otherwise identical to the transition with 1 elec- 
tron per atom, viz., bond dimcrization occurs along the 
diagonal directions, and OO and bond distortion rein- 
force each other cooperatively. 



B. | electron per atom 

In CuIr 2 S4 and MgTi 2 04, the distortion along the 
chain directions is not bond dimcrization but a period 
4 distortion 3 -. This is as expected for a Peierls transition 
in a ID chain with carrier density 0.5. We have therefore 
performed self-consistent calculations within our model 
Hamiltonian also for density 0.5. 

Not surprisingly, we do obtain self-consistently an or- 
bitally ordered state here for nonzero g with a = 0. Even 
within an essentially 100% orbitally-ordered state (mi- 
nority orbital charge density < 0.01), however, and with 
very strong a, we were unable to obtain a stable Peierls- 
distorted state. In all cases our self-consistent simula- 
tions converged to states with disordered bond distor- 
tions and charge densities, indicating vanishing bond- 
charge distortion in the thermodynamic limit. We ob- 



IV. DISCUSSION 

In summary, we have carried out numerical studies on 
the 2D checkerboard lattice with two degenerate direc- 
tional orbitals per site - a model system that like the 
spinel compounds can in principle exhibit OO-driven 
Peierls bond distortions and charge ordering in multi- 
ple directions. In addition to OO and bond modulation 
terms, our model Hamiltonian includes both intra- and 
inter-orbital e-e correlations that were treated within the 
UHF approximation. Although some of our results have 
strong implications for the spinels, it is useful to pre- 
cisely understand the differences between the spinel and 
checkerboard lattices such that the applicability as well 
as limitations of our model can both be understood. One 
difference between the two lattices is that the plaquettes 
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in the 2D checkerboard lattice do not correspond to the 
tetrahedra in the spinel lattice because of the difference 
between horizontal and vertical bonds in Fig.QJa) on the 
one hand and the diagonal bonds on the other — . What 
is more important in the present context is that the 00 
in our model is not driven by a band JT transition that 
destroys the degeneracies of the atomic orbitals in the 
model of [l3|. Within our model the two orbitals of dif- 
ferent symmetries on a given atom are both potentially 
active orbitals. 

We believe that our demonstration of the co-operative 
interaction between 00 and the Peierls instability in Sec- 
tion IIII Al where each broken symmetry enhances the 
other, is of direct relevance to the t2 9 -based spinel sys- 
tems, where qualitative discussions have suggested simi- 
lar results^ 3 -. Similarly, our observation that the period 
2 bond distortion persists for 1 electron per atomic site 
even for incomplete 00, with majority charge density 
as low as 0.8 per orbital, may also be of significance for 
the spinels. This should be particularly true for nonzero 
onsite Hubbard interaction, which will tend to decrease 
the amplitude of the 00. Complete 00 in the real sys- 
tems CuIr2S4 and MgTi204 requires that the energy gap 
due to the JT distortion is significantly larger than the 
Hubbard interaction. It is at least equally likely that the 
commensurate charge and bond distortions found in the 
experimental systems are not due to complete 00 but 
are consequences of the complex interactions between the 
crisscrossing chains in the spinel lattice, as in the checker- 
board lattice. 

The implication of the absence of the Peierls instability 
for the case of \ an electron per atom in the checkerboard 
lattice within one-electron theory is less clear. One possi- 
ble implication is that our results for the 2D checkerboard 
lattice are irrelevant for the three-dimensional (3D) spinel 
lattice because of the fundamental difference between 
them that has already been pointed out in the above. In 
CuIr2S4, the only active orbital following the 00 is the 
d xy orbital 13 , which has been excluded within our model. 
It is thus conceivable that the ID character of the active 
orbitals in CuIr2S4 following 00 is much stronger than 
in the checkerboard lattice, and this is what drives the 
metal-insulator transition in the real system. It is, how- 
ever, equally likely that the 3D interactions between the 
djcy-based chains are as strong as the 2D interactions in 
the checkerboard lattice (recall, for example, that com- 
mensurate periodicity for independent ID chains requires 
complete 00, see above). In this case our null result for 
the uncorrelated checkerboard lattice would imply non- 



negligible contribution of e-e interaction to the metal- 
insulator transitions in CuIr2S4 and LiRli204. Further 
theoretical work based on the 3D pyrochlore lattice as 
well as experimental work that determines the extent of 
00 in the real systems will both be necessary to com- 
pletely clarify this issue. 

Finally, assuming that e-e interactions play a role, 
which is subject to further investigations as discussed 
above, this raises an interesting question, viz., what ul- 
timately is the driving force behind the metal-insulator 
transitions in CuIr2S4 and LiRIi204? Three of us have 
argued elsewhere that for carrier concentration precisely 
0.5, there is a strong tendency to form a paired-electron 
crystal (PEC), in which there occur pairs of spin-singlet 
bonded sites separated by pairs of vacancies^. This ten- 
dency to spin pairing is driven by nearest neighbor an- 
tiferromagnetic (AFM) correlations (as would exist in a 
large Hubbard-U system) and is enhanced in the pres- 
ence of lattice frustration. Although the original calcu- 
lations arc for the anisotropic triangular lattice with a 
single orbital per site, the same tendency to spin-singlet 
formation can persist also in the spinel lattice. If the 
insulating state in the spinels CUIX2S4 and LiRli2 04 can 
be understood as a PEC with spin-singlet pairing driven 
by AFM correlations, it may further indicate that e-e 
interactions play an important role in superconductivity 
found in several structurally-related spinels. Whether or 
not e-e interactions play a role in the observed supercon- 
ductivity in the spinels LiTi204, CUPJ12S4, CuRli2Se4 
has remained a lingering questio n 32 ' 33 . If the insulating 
state in this class of materials is indeed a PEC, the super- 
conducting spinels should perhaps be included among the 
systems in which superconductivity is driven not entirely 
by BCS electron-phonon coupling. 
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